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Abstract. We investigate the magnetic response in the quantum spin Hall phase of the layered 
Kane-Mele model with Hubbard interaction, and argue a condition to obtain the Meissner effect. 
The effect of Rashba spin orbit coupling is also discussed. 



The Kane-Mele (KM) model is proposed to describe the quantum spin Hall effect (QSHE) of 
electrons on the two-dimensional honeycomb latticed] |2]. Here, we will show that, in a certain 
parameter region, the Meissner effect is obtained from the effective field theory of the layered 
KM model with an electronic correlation [3]. We use h = c = 1 unit and the Minkovskian metric 
gfj-u _ (ii(ig{i^ —^)-, where fi,u = 0, x, y. 

One of the essential ingredients of the KM model[Tl[2] is the intrinsic spin-orbit (SO) coupling 
Asoi which gives an effective magnetic field depending on spin and also an excitation gap 
A = 3\/3Aso to the electron[2]. Thus, the term yields quantization of the spin Hall conductivity 
(SHC)II1[2], 

where d is the interlayer distance. The model can also have the Rashba extrinsic SO coupling 
Xr, which breaks the inversion symmetry and is induced by an electric field perpendicular to 
the honeycomb lattice plane. The term also breaks the conservation of electron spin We 
consider Xji = case, first. 

We add the on-site Coulomb repulsion U. The Hamiltonian per layer is, 

H = c-Cj + ^Aso J2 ^iA^zCj 

<ij> 

+UY,nitnii, (2) 

i 

where Cj (c|) is the annihilation (creation) operator of an electron with spin at the i-th site and t 
is the nearest neighbor hopping. The second term is the intrinsic SO term consisting of the next 
nearest neighbor hopping, and Uij = -^{di x dl2)z = ±1, where di and d2 are unit vectors along 
the two bonds where the electron moving from site j to i passes. We assume that U is lower 
than the critical value Uc, above which the system tends toward the other insulating phase. [3] 



Namely, we consider the phase which is continuously connected to the quantum spin Hall state 
U = 0. 

Let us discuss how to deal with the electron correlation U. On-site coulomb repulsion can be 
written by the on-site spin-spin interaction; 

Uui^rii^ = ^{ni-f + nil) " -^{4^Cif- (3) 

The first term just gives the renormalization for the chemical potential and can be neglected. 
We introduce the auxiliary field ipi, which is a three-component vector in the spin space, and 
use the Stratonovich-Hubbard transformation [3] 



where, 



H ^ HsH = H + AH, 



= JY.i4sc^-^0^)^ (4) 



and therefore, 



HsH = t^4(^j + ^^so 

<ij> 

-II'^i-4|ci + ^^|(^ip. (5) 

i i 

The spin-spin interaction is eliminated in appearance, but we have a coupling between cpi and 
the electron spin, and a quadratic term of (pi instead. 

We consider the continuum limit and take into account the low energy electronic excitations 
around K and K' points in the Brillouin Zone [U [2]. We introduce the electromagnetic U{1) 
gauge field [5] and SU (2) spin gauge field via covariant derivative, 

iD^ = idfj, - eA^ + ■ -, (6) 

where clq = (p (the auxiliary field in the continuum limit) and a is a constant external field 
introduced artificially to estimate the spin current. We define a parameter 

AUa^d 

9 = > 0, (7) 

where a is the lattice constant, and the microscopic Lagrangian density is [6], 
£ = |^2)o - iv{DxTz(Jx + Dy(7y) + Ar^fj^Sx} ^ 

+ ^l«oP, (8) 

2 2/io 2g 

where ^ = ^ras is the eight-component Fermion field which is labeled by the eigenvalues of 
the diagonal components of valley spin r, sublattice spin a and real spin s/2. eg and denote 
the dielectric constant and magnetic permeability, respectively. Note that, except for the last 
term, the Lagrangian ([8]) possesses the ?7(l)em x U{\)z local gauge symmetry. The SU{2) gauge 
symmetry is broken down to U{l)z, since the SO term contains Sz- 



The calculation which will be shown below is quite similar to one presented in Ref.[7], although 
the physical meaning of the spin gauge field is different. We integrate out ^ and obtain the 
one-loop effective Lagrangian for the gauge fields in the low energy and long wavelength region. 
The result is [7], 

-^cflF = — -^-flo^ + Andi (9) 

And = ^ye^'^aldpA^+'-^-^ + ^-^iValf 

+(terms independent of and A^), (10) 

where Cmd stands for the induced part of the effective Lagrangian. The first term in Eq. (jlOp 
is the BF term[8]. The coefficient is the quantized SHC given in Eq. ([1]). Note that only 
couples to the electromagnetic gauge fields. This comes from the fact that the SU (2) symmetry 
is broken down to U{l)z symmetry by the SO coupling. Maxwell term is renormalized as[9l[7j 

e = eQ + 5e 

i - 1 I e'"^ ,12) 

/i /xo 67r|A|(i 

The elastic term for ao(= f/?^) is also induced. We can recognize that any potential terms (i.e., 
zeroth-order terms with respect to the derivative d^) of and also in C\nd are absent because 
ofC/(l)em X ^^(1)2 gauge symmetry in the Fermionic part of microscopic Lagrangian Thus, 
the low energy and long wavelength physics of and is described definitely by Eq. ([9]). 

The solution of the static equation of motion obtained from the Lagrangian ([9|) has been 
investigated |10^ [TT] . The argument indicates that when the dimensionless parameter satisfies 

<lt^g > 1, (13) 
we obtain the perfect diamagnetism as the energetically favorable solution. Namely, we obtain 

B{x) = B(0)-^e-'^o^, (14) 

where, 

1 



^0 ^ 



2ea'^y ^ fi afyfig - 1 ' 
shows the penetration depth. We also obtain 

«^(^) = 4^e--o^ (16) 

We can see from the BF term in eq. (|1U|) that electric current fiows perpendicular to the gradient 
of ao{x) jllj . and this current screens the applied magnetic field. 
Let us take into account the Rashba coupling. [H [2] We add 



(17) 



to the microscopic Lagrangian ([8]). Since the Rashba term breaks U{l)z symmetry, the 
symmetry-breaking terms are induced additionahy in the effective Lagrangian Eq. Up 
to the Gaussian approximation, the important alteration is the quadratic term of Oq; 



,22 ,1 I ^2 1 \ „22 



where [7] 



Eq. (|18|) indicates that the parameter g is renormalized as 

9 ^ . ^ 2 - (20) 

The other parameters a^y, e and in Eq. Q also receive renormalization of the order 0(A|j)[7j- 
The essential point is that Eq. (jl3p . which is the condition for the Meissner effect, is changed 
as (see eq. (pOj) ) 

> 1. (21) 



1 -5(771^ 



Compared with the condition (|13p . Eq. ()2ip becomes more attainable. gvrP'^ is positive and 
proportional to A^, which can be changed by the electric field perpendicular to the two- 
dimensional plane. 
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